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Abstract 



We present results for the bb spectrum obtained using an 0(MbV 6 )-correct 
non-relativistic lattice QCD action, where Mj, denotes the bare b-quark 
mass and v 2 is the mean squared quark velocity. Propagators are evalu- 
^ ' ated on SESAM's three sets of dynamical gauge configurations generated 

with two flavours of Wilson fermions at j3 = 5.6. These results, the first 
of their kind obtained with dynamical Wilson fermions, are compared to a 
quenched analysis at equivalent lattice spacing, [5 = 6.0. Using our three 
sea-quark values we perform the "chiral" extrapolation to m e g = m s /3, 
where m s denotes the strange quark mass. The light quark mass depen- 
dence is found to be small in relation to the statistical errors. Comparing 
the full QCD result to our quenched simulation we find better agreement 
of our dynamical data with experimental results in the spin-independent 
sector but observe no unquenching effects in hyperfine-splittings. To pin 
down the systematic errors we have also compared quenched results in 
different "tadpole" schemes as well as using a lower order action. We find 
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that spin-splittings with an 0{M\ ) v A ) action are O(10 %) higher compared 
to C(MfeW 6 ) results. Relative to the results obtained with the plaquette 
method the Landau gauge mean link tadpole scheme raises the spin split- 
tings by about the same margin so that our two improvements are opposite 
in effect. 
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I. INTRODUCTION 



Not just since the formulation of Heavy Quark Effective Theory a few years ago 
0, systems containing one or two heavy quarks have been the focus of much attention, 
both theoretically and experimentally. Today, B mesons are at the center of experimental 
efforts to determine quantities related directly to the CKM matrix, always in the hope 
that new physics might emerge. But, just as for systems with light quarks, the properties 
of hadrons containing heavy quarks are determined by non-perturbative physics and here 
the lattice can provide the missing link. 

Non-relativistic QCD (NRQCD) || is an effective theory used to describe heavy 
quarks, M q ^> Aqcd- The NRQCD Lagrangian is written as a series of opera- 
tors expanded in powers of the mean squared heavy quark velocity v 2 , £nrqcd = 
J2i,n c i(9 2 )Oi(M g v 2n ) where the coefficients q are obtained by perturbative matching 
with QCD. The lattice version of NRQCD |3|,[|] allows one to simulate heavy quarks 
with lattice spacing errors of 0(pa), where p is a typical momentum m Aqcd 5 rather 
than the usual 0(M q a) for Wilson-type heavy quarks. This makes lattice NRQCD a 
promising technique to simulate systems containing a 6-quark. 

Lattice NRQCD has been widely used in the past few years to calculate a variety of 
phenomenologically interesting quantities (see for a recent summary). The first step 
in the NRQCD program, naturally, is the calculation of the spectrum of the bb system. 
Such calculations have been performed using the quenched approximation by So 
far, there has only been one analysis of the T with dynamical quarks: the authors of 
PJ9|JT0|] applied an 0(M b v 4 ) correct NRQCD action to gauge fields incorporating the ef- 
fects of two dynamical staggered fermionsQ. It was found that the experimentally known 
spin-independent spectrum could not be reproduced using the quenched approximation 
whereas the data with 2 dynamical fermions agrees much better. These measurements 
were pushed to high statistics (4% statistical errors) in PJTDfl and used to obtain preci- 
sion measurements of the strong coupling constant a s . Quenching, the approximation in 
which the effects of dynamical fermions are disregarded, was concluded to be the largest 
source of uncertainty in this calculation. The dependence of the splittings on the dy- 
namical quark mass was estimated, using just two mass values, to be of the same order 
as the statistical error. No simulation exists, so far, applying NRQCD to dynamical 
Wilson fermions. 

In the spin- dependent sector the emerging picture is much less clear. The P-fine 
structure in the unquenched theory [j5|,|TQH , obtained with the (9 (M;,?; 4 ) -correct action, was 



found to be in very good agreement with experiment whereas the quenched results, using 
the same action, predict much smaller splittings. However, recent results of ||ll,12j] (all 



of which use the quenched approximation) have exposed the sensitivity of the spin- 



lr rhe high statistics measurements in the spin- independent sector in were obtained using 
an 0{Mi J v 2 ) correct action. 
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dependent splittings to the details of the action: several improvements, such as the 
inclusion of higher order spin-dependent terms in the NRQCD expansion, the addition of 



lattice-artefacts correcting terms of Oia ) ||8|JT^1 as well as an improved phenomenological 



estimate of the coefficients q [|12|)[yl were found to have sizeable effects of O(10 — 20%) 
for these splittings. Thus, there has recently been strong support in favour of using the 
0(Mf,f 6 )-correct NRQCD action, where, at highest order, all spin- dependent corrections 
are added, as well as for the use of the mean link Landau gauge tadpole scheme to 
estimate the perturbative coefficents c,. The effects of these improvements for dynamical 
configurations have not been studied. 

In the simulation presented here we attempt to address many of the open issues pin- 
pointed above. Using SESAM's large sample of dynamical Wilson-fermion gauge con- 
figurations we study both the spin-independent as well as the spin- dependent spectrum 
of the T. Our strategy, in searching for sea-quark effects, will be to compare our final 
dynamical results to that of a quenched simulation at equivalent lattice spacing. Thus, 
we hope to consolidate that unquenching brings the spin-independent splittings into 
good agreement with experiment. Using our three sea-quark masses we can also study 
the dependence of mass-splittings on the light sea-quark mass. Following the recent 
suggestions of H8|jl2]| we have implemented the NRQCD action including spin-dependent 
corrections of 0(Mf,v 6 ) and we remove tadpoles using the mean link calculated in Lan- 
dau gauge. With these ingredients we hope to clarify the effect of unquenching in the 
spin-dependent splittings. In addition, using the quenched configurations we investigate: 
(i) the effect of changing tadpole prescriptions; (ii) the effect of using an 0(M b v*) correct 
action compared to the 0(Mf)V 6 ) corrected one. 

II. SIMULATION DETAILS 

In section [11 A| we present the discretized version of our (9(M{,t> 6 ) correct NRQCD 



action. Section [II B| addresses issues concerning the determination of the perturbative 
coefficients Ci(g 2 ). 



A. Action 

The non-relativistic (Euclidean) lattice Hamiltonian to 0(MbV 6 ) consists of ||: the 
kinetic energy operator, 

A( 2 ) 

H = -=— , (1) 
which is of order M^v 2 ; relativistic corrections of order M^w 4 , 



SH^ = -ci^— f- + c 2 -^ (A -E-E- A) 



3M 6 3 8M b 2 
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-cs-^ra • (A x E - E x A) - c 4 -?—a ■ B 



a 2 A( 4 ) a(A (2) ) " 
24M ft C6 16nM| 



and spin-sensitive corrections of order M^v 



^ (6) = -^{A (2) >-B} 



M3 



-c 8 -^-{A^, cr • (A x E — E x A)} 



(2) 

64M b 4 1 

. 2 

-cg-^TT^ ■ E x E . (3) 

We do not include spin-independent corrections of 0(MbV 6 ). The E and B fields are 
definded in terms of the field strength tensor 

F^(x) =--^ I -Tr(lm[/^(x)) I , 
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5* = --e^F^ , (4) 

where, is the standard plaquette and the sum is over all anticlockwise plaquettes in 
the /x-i/-plane, 1 < /i < 4, \ < u < ji. 

The removal of (9 (a 2 MbV 4 ) discretization errors is achieved by: (i) using an improved 
version of the lattice field strength tensor || in the leading spin-dependent interactions 
(c 3 and c 4 ): 



3 p 6 

+£/J(x - fi)F^(x - - A) - (/x «-> i/)) ; (5) 

(ii) adding the last two terms in <5if ( - 4 ' 1 to (classically) correct for finite lattice spacing 
errors in the spatial and temporal derivatives and (iii) by using the improved derivative 
operator 

K = a, - |a«aa h (6) 

aA^ ] G(x) = Ui(x)G(x + at) - G(x) 
aA\~ } G{x) = G{x) - U}(x - ai)G(x - ai) 

aAjG(x) = - (jJi(x)G(x + ai) — Uj(x — ai)G(x — ai)^j 



5 



in the leading spin-dependent interactions (c 3 ). 

For completeness we also give our definitions of A^ 2 ^ and A^: 

A ( 2 ) =EA p) 

i 

a 2 A^ } G(x) = Ui(x)G(x + at) + U}{x - ai)G(x - at) - 2G(x) 

i 

Following |7]] the quark Greens functions are calculated from the evolution equation 

aH \ n x ( ^ aH \ n 



G(t+l)=(l--^j U l[}-^) (l-aSB)G(t), (8) 
G(0) = 5 x ,o , 

where G(t) is a two-spinor component object. 

The parameter n serves to stabilise the evolution for small bare quark masses and is 
set to n = 2 01. 



B. Tadpoles 

An important feature of our simulation is the use of "tadpole" -improvement: gauge 
links Up are replaced by U^/u where the mean link u accounts for the effects of tadpole 
diagrams |TJ]JTH]. The coefficients q are then set to their tree-level values of 1. 



The choice of uo is, of course, not unique; here we choose uo to be the mean value of 
the link variable in the Landau gauge: 

4 L) = (\™,), d,A, = 0. (9) 

With this choice of gauge (TrC/^) is maximised so that u is as close to 1 as possible; 
it was therefore argued in JT3] that remaining tadpole contributions cannot be due to 
lattice artefacts. A recent NRQCD lattice calculation of the cc spectrum over a variety 
of lattice spacings suggests that this choice improves the scaling behaviour of spectro- 



scopic quantities fll2[ , notably the hyperfme splittings. This view is also endorsed in 
PJTT[| where the authors have studied the scaling behaviour in the quenched bb spectrum 
using three values of the coupling (/3 = 5.7, 6.0, 6.2), an 0(M b v 4 )-correct action and the 
plaquette tadpole value (defined below). They find the spin- independent spectrum to 
display insignificant scaling violations whereas the hyperfme splittings do not scale that 
well. This is attributed partly to 0(a 2 ) errors in the B fielclQ arising in the term a-B but 
partly also to the choice of plaquette tadpole improvement not capturing the tadpole 



2 which was obtained fom the unimproved field strength tensor, eq.|j. 
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effects sufficiently well. Additionally, two more arguments in favour of are given 
in fl3f : firstly, the static potential shows less violation of rotational invariance using 



a w^-tadpole improved gluonic action and, secondly, the non-perturbative determina- 
tion of the clover coefficient [16|] is in good agreement with the w^-tadpole improved 
perturbative result. Since, for the bb, all simulations performed to date, notably those 
with dynamical staggered quarks, have been using Uq = (^TrU^)^, our measurement 
will help to explore the size of the systematic error associated with the choice of u . 
In the quenched sector we have also performed an analysis of the S-splittings (where 
signals are cleanest) using the plaquette prescription and, furthermore, we compare the 
(9(Mf,f 6 )-results to a simulation where the NRQCD action is correct to 0{M\ ) v i ). In 
table | we compare our values of to the commonly used Uq P \ Note that uq weights 
the field strength tensor with four powers so that, naively, using Uq P ^ instead of Uq may 
change the hyperfine splittings by as much as 8%. 



C. Correlators and Smearing 



Meson correlation functions at zero momentum are built from quark propagators 
combined with suitable interpolating operators: 



G%S? (*, P = 0) = E Tr [Gt ( X , t ) T\ sk) (y - x) G (y, t) 



(10) 



where the source smeared propagator G is obtained by evolving the extended source: 

G(y) = £G(y-x,t)r (sc) (x) (n) 

X 

The interpolating operator T^ sk \x) = Q$( sc / sfc ) (x) contains a spin matrix f2 and a 
spatial smearing function $; the superscript sc/sk denotes source and sink smearing re- 
spectively. A proper choice of $ is crucial to obtain clear signals for excited states. Here, 
we benefit from a recent lattice determination of the bb potential including relativistic 
corrections ||17|| . These authors start from the 0(v 4 ) non-relativistic Schrodinger-Pauli 
Hamiltonian fIB |, which, in the center of mass frame of the quark and antiquark, has 
the form 



H 



i=l,2 V 



+ J^V 2 (F (r) + V a {r)) 



2rrii 8m'- Sm. 
+V Q (r) + V sd (r, L, S x , S 2 ) + V vd (r, p) 
It contains the central potential of Cornell type: 

V (r) = V + kr - - , 
r 



(12) 



(13) 



spin- and velocity-dependent (sd, vd) potentials (see [0) and the Darwin-term V a . 
Following W71 we numerically integrate the radial Schrodinger equation 
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9E, L (r) = F(r, E, L, a, m, k) g EtL {r) (14) 
for definite radial quantum number and angular momentum, where 

F(r, E, L, a, m, k) = L ^t -m(E-kr+-\ . (15) 

The term — - contains the Coulomb-term of the central potential, the - contribution of 
the Darwin term as well as the leading - contribution from Vvd- We leave m and k as free 
parameters but fix the contributions of the terms not coming from the central potential 
to the lattice values found in ]17|]. The Coulomb coefficient e is that of the quenched 



simulation in [I7| augmented by 20 % to take into account the different running at short 
distances observed in A direct comparison with hydrogen wave functions clearly 
demonstrates the superiority of our choice. 

The continuum wave functions are converted to the lattice using the scales obtained 
from SESAM's p-masses for each K sea . The spin matrices Q are taken from [0] and listed 
in table |TJ. Note that the finite difference operators included in table |H| f° r t ne P- waves 
are only applied in the case of delta-function smearing. 

For the T and r]b we calculate a 4 x 4 matrix of correlators with four different 
smearings at source and sink, sc/sk = local, 1,2, 3 corresponding to a point source 
(local), the ground state (1), the first (2) and second (3) excited states respectively. For 
the L = 1 states we restrict ourselves to the ground state and the first excitation as 
signals deteriorate. Correlators with momenta upto |p| = 2 are also calculated. Gauge 
configurations are fixed to Coulomb gauge. 

D. Lattice Parameters 

The lattice parameters we have used are displayed in table [[II]. SESAM has re- 
cently completed the generation of gauge configurations with two (degenerate) dynamical 
fermions at three values of the sea-quark mass. For each mass the run consists of (9(5000) 
trajectories^ and we analyse every 25th trajectory giving a total of 200 configurations at 
each sea-quark. This is by far the largest set of dynamical Wilson-fermion gauge config- 
urations to date, more than a factor two larger than the sample analysed with respect to 
NRQCD in HQ. Details concerning the generation of our dynamical configurations and 
issues surrounding autocorrelations are discussed in ref. ||20|| . Here, we mention only the 



values of ir to p mass ratio and the lattice spacing obtained from the p - which we have 



measured in a standard light spectrum calculation on the lattice fl2T|,|2"2"| - for our three 



3 In a standard Hybrid Monte Carlo algorithm we typically set the time-step size dt = 0.01 
and perform 100 molecular dynamics updates to generate one trajectory. 

4 The authors of [||] analyse 400 configurations taken from a total trajectory length of 2000. 
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dynamical sea quarks K sea = {0.156,0.157,0.1575}: ^ = {0.839(4), 0.755(7), 0.69(1)} 
and a- 1 = 2.33(6)GeV. 

Taking advantage of the smallness of the bottomonium system, we exploit configu- 
rations more than once by starting the propagator evolution both at different spatial 
source points and on different timeslices. A binning procedure confirms that our 4 
measurements per configuration are indeed independent. 

Errors are obtained from correlated fits using the bootstrap procedure. They corre- 
spond to 68 % confidence limits of the distribution obtained from 300 bootstrap samples. 

Throughout, we fix the bare heavy quark mass to a value aM& = 1.7. In a quenched 
simulation by the NRQCD collaboration at (3 — 6.0, a bare b quark mass of 1.71 was 



found to yield the correct physical value of the T [23 



III. SPECTROSCOPY 

We adopt spectroscopic notation and denote (radially excited) spin-parity eigenstates 
by n 2S+1 L j (see table 0). 

Figure [3] shows a subset of local masses: although we have not tuned the smearing 
functions we generally find good plateaus for ground and excited states. For L = the 
radially excited correlators remain in the first/second excitation for about ten timesteps. 
Local masses for P states are noisier and drop to the ground state without staying in an 
'excited plateau' as nicely as for the S states. 



A. Fitting procedures 

To extract energies we fit several smearing combinations of correlators simultaneously 
to a multi-exponential ansatz. We find that vector fits to smeared-local propagators, 



G^ n (t) = £ b SCjk e- E ^ , (16) 

k=l 



are quite stable whereas matrix fits appear to require higher statistics. Energy levels are 
obtained from correlated two-exponential fits to the two smearing functions sc = 1,2 
with eq. |l(| Hyperfine splittings are determined from single exponential fits to the ratio 

^mesoriA 

( sc > sk ) (f\ - A p -(E A -E B )-t / 17 a 

^mesoiiB \ L J ~ ^ c ; \ L 1 ) 
Lr (sc,sA:) 

thus exploiting the strong correlation between the data. We use smeared-local and 



smeared-smeared correlators in eq. [L7]. The results quoted correspond to the value with 
the lowest Q-fit-parameter |24| but with the statistical error enlarged to take into account 
the uncertainty in the fitting ranges, thus encompassing several possible fit-results from 
different fit-ranges. For the hyperfine splittings the errors allow for the spread in fit- 
results obtained using smeared-local as well as smeared-smeared data. More complicated 
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fit-functions confirm the results obtained from the single exponential ansatz but do not 
behave as stable. Results of our fits are given in tables [TV| and |V} 



B. Light quark mass dependence 

We now turn to the dependence of splittings on the dynamical quark mass. Following 
the analysis in ref. [25J we extrapolate energy level splittings linearly in the quark mass: 



Am = Amo + c 



m 



u,d,s 



where m q denotes the bare light quark masses. The relevant scale we pick is given by 
m u +m d +m a ^ corresponding to a "physical pion mass" of m ^ + ^ m i< |25|]26|| . Since gluon 
momenta inside the bb are much larger than any of the three lightest quark masses, T- 
physics should be insensitive to the m u 4 — m s mass-splitting. The value of the strange 



quark mass is taken from our recent light quark mass calculation |2T| and is such that 



m cS = ^ = ^(---) ^0.0156 (19) 

6 OCL \K. K ( 



lies close to our lightest sea-quark mass, am q (K = 0.1575)0. The authors of [2S|] have 
noted that extrapolating to m cS = ^ is appropriate for the 3-flavour case. In this 
context, we point out that there is an uncertainty in the strange quark mass m s due to 
the fact that we have simulated only two (degenerate) dynamical flavours. These are 
identified with the up and down quarks so that, effectively, m s is defined in the 2-flavour 
theory in a "partially quenched" way only pl| . Thus, it is not clear what precise value 
of m e s mass splittings should be extrapolated to. 

The extrapolations are shown for various splittings in figures |2| and [j], table |V] 



contains Am(m s /3) and Am . No significant deviation from the linear parametrization 
is apparent in the data. With the exception of the 3 Si — 1 S splitting the light quark 
dependence is too small to be resolved given our current sensitivity. Within large errors, 
we find Am l2 s/3) , which can be taken as an upper limit of the uncertainty in m e s, to 
be 3 — 20% for the spin- independent splittings. Choosing m e s < m s /3 would have the 
effect of increasing the inverse lattice spacing a -1 (see next section). Clearly, much 
higher statistics is necessary to pin down these effects. 



C. Lattice spacing 

The lattice scale is taken from the average of the 2 3 S'i — l 3 Si and CG( 3 P) — l 3 Si 
splittings at m s /3 (table [VI I|) . Recall that these splittings are virtually insensitive to 



5 The three sea-quark values correspond to pions am n = {0.445(2), 0.341(3), 0.276(5)} J22j| . 
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the heavy-quark mass and vary only by a few percent from the bb to the ccf]. Note that 
the average lattice spacing for j3 = 5.6 for n/ = 2 agrees very well with the quenched 
one at /3 — 6.0 so that we can compare our results in the "full" and in the quenched 
theories. 

The ratio of splittings, R SP = (2 3 S 1 - l 3 S 1 )/(CG( 3 P) - l 3 ^) is shown in figure |(a) 
as a function of the light quark mass. All values agree with the experimental number of 
1.28. Finally, we use the average lattice spacings - at m s /3 in the full theory - to convert 
the lattice numbers to their physical values, see table |V111| . 



IV. DISCUSSION 



Figure [| summarizes our results for 2 and flavours from table |VIHj Unquenching 



effects are clearly visible in the spin-independent part of the bb spectrum; this is also 
evident from table [VI]| which shows that the lattice spacings from the IS — 2S and the 
IS — IP splittings do not agree in the quenched theory whereas they coincide when 2 
dynamical quarks are switched on. The ratio of these splittings is plotted in figure |](b) 
as a function of the number of flavours. 

We do not observe any significant impact of unquenching on the hyperfine splittings. 
In particular, the P-hyperfine splittings seem to be underestimated for both rif = and 
n/ = 2. Clearly, this result needs to be corroborated with higher statistics. Errors on 
our hyperfine splittings encompass the statistical error as well as the uncertainty in the 
fitting range; the latter is essentially of the same size as the statistical error so that we 
can expect some improvement with higher statistics. 

Compared to the only other simulation of NRQCD with dynamical quarks, performed 
by the NRQCD collaboration^ using staggered fermions [0,0], we have introduced two 
new features. Firstly, we have included all spin-dependent corrections to 0(M^v & ) and 
secondly, we have chosen an alternative way of removing the tadpole diagrams^. Com- 
pared to an 0{Mi)V A ) correct action we expect the inclusion of 0(Mf,v 6 ) terms to effect 
the hyperfine splittings on the 10 % level, since v 2 ~ 0.1. Performing a relatively cheap 
0(M h v A ) quenched simulation at (3 = 6.0 (with 200 measurements) and calculating only 
the 3 Si — 1 Sq and the (spin-independent) 25 — IS splittings we find this naive expecta- 
tion to be very well satisfied as is shown in figure || Very recently, a similar effect was 
reported in || (0(10 %) for the bb) as well as in fll2"] , where the effect of next-to-leading 



6 The 2 3 Si-l 3 Si splittings are 562.9 and 589.1 for the bb and cc respectively; for the CG( 3 P)- 
l 3 Si they are 439.8 (bb) and 428.4 (cc) - all in MeV j27|. 



7 on gauge configurations generated by the HEMCGC group. 

8 Other minor differences in the two simulations are: (i) a bare 6-quark mass of 1.8 was used 
by the NRQCD collaboration; (ii) the HEMCGC data was not "chirally" extrapolated. 
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order spin-dependent interactions for charmonium was found to be of O(60 %). Ref. || 
also found a rather severe shift, due to the 0{M\,v 6 ) terms, of O(30 %) in the P-fine 
structure of the bb. These results confirm Trottier's conclusion [|l!| that the first three 
terms in the expansion for the hyperfine splitting are oscillating in sign. Therefore, the 
use of an 0(Mi,v 6 )-correct action is mandatory to get meaningful results for the fine- 
structure of the T. 

In addition to the order in the NRQCD expansion, the hyperfine splittings will be sen- 
sitive to the value of u . Naively, one expects the splitting to vary as 1/uq (since the 
B fields contain four links) and the ratio of the Uq with the plaquette prescription to 
that with the mean link Landau gauge value is roughly 10 % (see table | for values of 
Mo in both schemes). We therefore also plot in figure || the 3 Si — 1 So splitting obtained 
with the 0(M b v 4 ) correct action but using the plaquette prescription for u (this re- 
sult is from the NRQCD collaboration^) . The plot confirms the naive expectation and 
shows that changing uq to the mean link Landau value shifts the hyperfine splitting in 
the opposite direction as adding the 0(Mi,v 6 ) spin-dependent corrections. For reference 
the corresponding value is also shown in the spin-independent 25 — IS splitting where, 
within the larger error, no effect is seen. 

With this in mind, we compare in figure [F| our {0(MbV 6 ), tijj , rif = 2} results to the 
{0(Mi,v 4 ), Uq P \ rif = 2} of the NRQCD collaboration ||. Both in the spin-independent 
and in the spin-dependent sectors the results are in good agreement although our split- 
tings tend to be slightly smaller. 

V. SUMMARY AND OUTLOOK 

We have presented the first calculation of the bb spectrum with dynamical Wilson 
fermionsQ. Our non-relativistic lattice action is correct to 0(MbV 4 ) for spin-independent 
operators and includes all spin-dependent corrections of 0(Mi,v 6 ). We rely on tadpole 
improvement choosing the mean link in Landau gauge as our improvement scheme. 

By use of appropriate smearing techniques we obtain clean signals for ground and 
excited states. 

We have studied the light quark mass and flavour dependence of the bb spectrum. All 
quantities were linearly extrapolated to m s /3. We are helped in these extrapolations by 
the fact that m s /3 is very close to our lightest sea-quark mass so that the extrapolated 
values are easily consistent within errors with the values at k = 0.1575. The lattice spac- 
ings were determined after "chiral extrapolation" with an error of 6 %, which accounts 
for the statistical and systematic (fitting) uncertainties. 



9 We have reproduced their value on 200 configurations but with a larger error. 
10 We note that the authors of |2£| are presently starting a similar analysis. 
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Comparing to the quenched calculation at similar lattice spacing we find a distinct 
unquenching effect in the spin-independent splittings which can be quantified by the 
ratio Rsp = determined to be 1.23(11) with 2 dynamical flavours and 1.43(9) 

for the quenched case; the experimental value is 1.28. 

We have estimated the systematic error arising due to two different choices of the 
tadpole subtraction factor u ; whereas spin-independent quantities remain unaffected 
by changes in uq the fine splittings increase when using the mean link Landau gauge 
Uq. Interestingly, we find this effect to be of the same magnitude but of opposite direc- 
tion as switching on the spin-dependent 0(MbV 6 ) corrections. These two effects thus 
combine to yield good agreement between our {0(MbV 6 ), Uq, rif = 2} results and the 
{0(MbV 4 ), Uq P \ nf = 2} results of the NRQCD collaboration with dynamical staggered 
fermions ||. 

Our results confirm that the 0(Mi,v 6 ) terms in the action are mandatory for calcu- 
lating the hyperfine splittings. As expected, the determination of the lattice spacing a -1 
(as used, for example, in the determination of a s [Kj) remains unaffected by changing 
tadpole schemes. 

Overall, NRQCD has proven to work very well in the spin-independent sector, in 
particular in giving a precise determination of the lattice spacing a. Given the weak 
dependence of these splittings on the heavy quark (experimentally) as well as on the light 
quark masses, as found here, it appears worthwhile to push the scale determinations to 
higher statistics using our configurations generated by T%L at k — 0.1575 on 24 3 x 40 
lattices (j3 = 5.6). This would enable a reliable extrapolation to 3 or 4 flavours. The 
dominant source of error in such a lattice scale determination is most likely due to the 
remaining lattice discretisation errors in the gauge configurations. 

Although higher statistics are highly desirable, it seems unlikely that progress in the 
spin-dependent sector will come from this approach alone. It is worrying to find the 
spin-dependent corrections of 0(M b v 6 ) as large as 10 % even for the bb, decreasing the 
fine-splittings relative to the lower order estimates. Even at this order the NRQCD 
expansion does not describe the hyperfine structure of the spectrum. A perturbative or 
non-perturbative calculation of the coefficients Cj is badly needed^. 
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11 Work in this direction was recently started by [29]. Perturbative calculations [30| indicate 
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FIGURES 





FIG. 1. Two examples of effective masses with several radial excitations. The r]b is at our 
lightest sea-quark, k = 0.1575 whereas the hb is at the heaviest, k = 0.156. 
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FIG. 2. Extrapolation of spin independent splittings to m s /3. The IP — l^Si splitting is 
plotted for comparison as it is one of the splittings used to determine the scale (it is shifted 
downward by 0.1 on the right plot). The triangular symbol denotes the extrapolated value. 
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TABLES 

TABLE I. Comparison of tadpole improvement schemes. 



uo k = 0.156 k = 0.157 k = 0.1575 quenched 

u ( P = (^TrU^u) 0.8688 0.8695 0.8697 0.8778 

u ( L) = (\TrUp) 0.8499 0.8519 0.8529 0.8608 



TABLE II. Operators and naming scheme. 



bb 


2S+l Lj{J PC ) 


Q 


Vb 


^0(0-+) 


1 


T 


3 Si(l-) 


°i 


h b 


^1(1—) 


A, 


Xb 


3 Po(0 ++ ) 




Xbi 


3 Pi(l++) 


AiCTj — AjCTi 


Xb 2 


3p 2 (2++) 


AjCJj — Aj(7j 


X 


CG( 3 P) ~P= 





TABLE III. Simulation details (n/ denotes the number of flavours). 





/3 dyn = 5.6, n f = 2, 16 3 x 


32 






^sca 


0.156 




0.1570 


0.1575 


number of configurations 


200 




200 


200 


measurements 


800 




800 


800 




Ajuenched = 6 - ' n .f = °' 163 


x 32 








number of configurations: 


200 






measurements: 800 
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TABLE IV. Spectrum results in lattice units for all three sea-quark masses and for the 
quenched simulation. 









nf = 2,/3 = 5.6 




n f = 0, (3 = 6.0 


Level 




k = 0.1560 


k = 0.1570 


k = 0.1575 




1% 




0.3493(8) 


0.3476(7) 


0.3446(7) 


0.3309(8) 


2 1 S 




0.635(20) 


0.613(15) 


0.588(18) 


0.574(12) 


l 3 Si 




0.3651(9) 


0.3621(9) 


0.3582(9) 


0.3438(8) 


2 3 Si 




0.635(21) 


0.610(15) 


0.584(17) 


0.586(13) 


3 3 Si 








0.75(5) 


0.81(2) 


l^i 




0.559(11) 


0.548(9) 


U.04i(iU J 


u.oiz^y J 


2 1 P 1 




0.80(4) 


0.74(2) 


0. / 1(2 J 


0. ( 2(2 ) 


TABLE V. Splittings at fixed sea-quark as well as 


quenched results (in 


lattice units). 








n f = 2, (3 = 5.6 




n f = 0, (3 = 6.0 


Splitting 


k = 0.1560 


k = 0.1570 


k = 0.1575 




2^0- 




0.285(20) 


0.266(9) 


0.243(7) 


0.243(8) 


2 3 5i- 


-l 3 5i 


0.269(25) 


0.248(15) 


0.226(17) 


0.242(13) 


3 3 5i- 


-l 3 ^ 






0.39(6) 


0.47(2) 


l^i- 


-l 3 5i 


0.194(10) 


0.186(9) 


0.183(10) 


0.168(9) 


2 1 P 1 - 


-I'Pi 


0.24(4) 


0.19(2) 


0.17(2) 


0.21(2) 


1P- 




0.194(11) 


0.186(9) 


0.182(10) 


0.169(9) 


l 3 5i - 


-l l S 


0.0160(2) 


0.0143(2) 


0.0138(2) 


0.0126(2) 


1 3 P 2 - 


-I'Pi 


0.0046(10) 


0.0031(10) 


0.0032(10) 


0.0050(6) 


l^i- 


- l 3 Pi 


0.0023(11) 


0.0028(9) 


0.0017(6) 


0.0016(8) 




-l 3 Po 


0.0105(19) 


0.0115(15) 


0.0115(11) 


0.0098(8) 


1 3 P 2 


- IP 


0.0039(8) 


0.0036(7) 


0.0033(5) 


0.0038(5) 


IP - 


1 3 P1 


0.0029(9) 


0.0023(7) 


0.0017(6) 


0.0027(7) 


IP - 


l 3 Po 


0.0111(15) 


0.0110(13) 


0.0115(11) 


0.0110(9) 


IP - 


1 X P1 


0.0006(7) 


-0.0005(7) 


-0.0001(5) 


0.0011(4) 
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TABLE VI. Splitting results (in lattice units) extrapolated to m s /3. 



Splitting 


Am(m s /3) 




2 iSn — 1 ^ iS'n 


0.238(9) 


0.212(18) 


2 3 S'i - l 3 5i 


0.223(18) 


0.201(29) 


l x Pi - l 3 Si 


0.181(9) 


0.176(16) 


2 1 Pi - l 1 Pi 


0.16(3) 


0.13(5) 


IP - l 3 5i 


0.181(9) 


0.175(16) 


l 3 5i - l 1 ^ 


0.0135(2) 


0.0124(3) 


1 3 P 2 - l 1 P 1 


0.0029(9) 


0.0021(15) 


1 X P! - 1 3 Pl 


0.0019(7) 


0.0015(13) 


1 x Pi - 1 3 P 


0.0117(13) 


0.0123(22) 


1 3 P 2 - IP 


0.0032(6) 


0.0029(9) 


IP - 1 3 P1 


0.0016(7) 


0.0010(11) 


IP - 1 3 P 


0.0114(12) 


0.0116(21) 


IP - 1 X P1 


-0.0003(5) 


-0.0007(10) 



TABLE VII. Determination of the lattice spacing from the 2 3 Si - l 3 5i and CG( 3 P) - l 3 Si 
splittings. We use the average value to convert our results to physical units. Psp is to be 
compared to the experimental value of 1.28. 





Splitting 


a-^GeV] 


Average a 1 [GeV] 


Rsp 


n f = 0, (3 = 6.0 


X-T 


2.33(10) 
2.61(13) 


2.47(10) 


1.43(9) 


k = 0.1560 


X-T 


2.09(18) 
2.27(12) 


2.18(12) 


1.39(13) 


k = 0.1570 


X-T 


2.27(15) 
2.36(10) 


2.32(10) 


1.33(12) 


k = 0.1575 


X-T 


2.49(17) 
2.41(10) 


2.45(12) 


1.24(13) 


m s /3 


X-T 


2.52(20) 
2.43(12) 


2.48(14) 


1.23(11) 
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TABLE VIII. Overview of our results in physical units for full QCD and the quenched case. 





Simulation Result [GeV] 


Experiment 


27] 


Splitting 


m s /3 


rif = U , p = o.U 




2 i^Vi — 1 i^Vi 
2 3 5 x - l 3 ^ 
3 3 5i - l 3 5i 

l 1 ^! - l 3 Si 
2 l P X - l 1 ^! 

IP - l 3 5i 


0.591(26) 
0.553(25) 

0.450(22) 

0.40(8) 
0.448(22) 


0.600(20) 
0.598(20) 

1.15(7) 
0.416(13) 

0.52(7) 
0.419(12) 


0.5629 
0.895 

0.4398 


l 3 5i - 1% 

1 3 P 2 - l 1 ^! 
l 1 ^! - l 3 Pi 
l 1 ^ - 1 3 P 


0.0334(19) 
0.0071(25) 
0.0046(18) 
0.0291(30) 


0.0313(12) 
0.0122(15) 
0.0039(20) 
0.0243(21) 




1 3 P 2 - IP 
IP - l 3 Pi 
IP - 1 3 P 
IP - l 1 ^! 


0.0080(15) 
0.0040(18) 
0.0283(35) 
-0.0009(15) 


0.0094(12) 
0.0067(15) 
0.0271(20) 
0.0028(10) 


0.0130 
0.0083 
0.040 
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